The Radiative Return: a review of experimental results 



Achim Denig a 

a Universitat Karlsruhe, Institut fur Experimentelle Kernphysik 
Postfach 3640, D-76021 Karlsruhe, Germany 

The radiative return is a new method for hadronic cross section measurements at electron-positron colliders, 
which are operated at a fixed center-of-mass energy (so-called particle factories). In order to lower the effective 
hadronic mass Mhadr only such events are taken, in which one of the electrons or positrons has emitted an initial 
state radiation photon. We present precision measurements of the pion formfactor from the Frascati (^-factory 
DA$NE with the KLOE experiment and measurements of higher particle multiplicities as well as a measurement 
of the timelike proton-antiproton formfactor from the BaBar experiment at the B-factory PEP-II. These raditative 
return measurements are compared to results, which are obtained by means of an energy scan, i.e. by means of 
a systematic variation of the beam energy of the collider. We also report on the impact of these measurements 
on the hadronic contribution of the anomalous magnetic moment of the muon, which is obtained via a dispersion 
integral using hadronic cross section data as input. 



1. THE RADIATIVE RETURN AND 
ITS CONNECTION TO THE MUON 
ANOMALY 

1.1. The radiative return method 

Modern particle factories, such as the Frascati 
0-factory DA$NE or the B-factorics PEP-II and 
KEK-B are designed for a fixed center-of-mass 
energy yfs. An energy scan for the measurement 
of hadronic cross sections is therefore not feasi- 
ble. A new and complementary ansatz is given 
by the 'radiative return' method, which utilizes 
events, in which one of the electrons (positrons) 
has emitted a photon. The invariant mass Mhadr 
of the hadronic system is reduced in this way and 
the hadronic cross section in the energy range 
Mhadr < y/s becomes accessible. In order to de- 
duce from the measured radiative cross section 
<r(e + e~ — + hadrons + 7) the non-radiative cross 
section er(e + e~ — > hadrons), a theoretical radia- 
tor function is used via the relation: 

, r 2 ^ cr hadrons+7 _ HY »/f 2 ^ 
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The radiator function, which describes the initial 
state radiation (ISR) process, has been com- 
puted up to NLO by the Monte-Carlo generator 
PHOKHARA [T] with a precision of 0.5%. If the 



measured radiative cross section is divided by the 
radiative muon cross section e + e~ — * /i + /i~7, the 
ratio R — er(e + e~ — * hadrons)/cr(e + e~ — * /i + 
can be deduced directly. In this case one does 
not need to normalize the data to the integrated 
luminosity, nor is the radiator function needed. 
However, an excellent separation of hadrons from 
muons is required. 

In the following we present results from the 
KLOE experiment at the Frascati collider 
DA$NE, which is operated on the 0(1020)- 
resonance and allows to measure the pion form- 
factor below 1 GeV. As we will see in the next 
subchapter, this measurement is of crucial im- 
portance for the hadronic contribution to the 
anomalous magnetic moment of the muon and a 
precision ~ 1% or better is needed. After that 
we discuss a series of results from the BaBar- 
experiment for the timelike proton-antiproton 
formfactor and for several exclusive final states 
with higher multiplicities in the mass range from 
threshold up to 4.5 GeV. The combination of 
KLOE and BaBar data allows to cover the 
hadronic cross section in the entire mass range 
below ~ 4.5 GeV. This is the relevant energy 
region needed for a significantly improved deter- 
mination of the anomalous magnetic moment of 
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Figure 1 . The contribution of different energy ranges 
of hadronic cross section data to the absolute value 
and the squared error for a^ dl and Aa^ dr (M|). 



the muon. 

1.2. The anomalous magnetic moment of 
the muon 

The recent measurement of the muon anomaly 
a M at the Brookhaven National Laboratory with a 
precision of 0.5 ppm [2] has led to renewed interest 
in accurate measurements of the cross section for 
e + e~ annihilation into hadrons. Hadronic contri- 
butions to the photon spectral functions due to 
quark loops are not calculable in the framework 
of perturbative QCD. It is well known, however, 
that the hadronic piece of the spectral function 
is connected by unitarity to the cross section for 
e + e~ — > hadrons. A dispersion relation can thus 
be derived, giving the contribution to as an 
integral over the hadronic cross section. 
Also the hadronic contribution to the running of 
the electromagnetic fine structure constant at the 
Z-pole Aa| i ^ ) dr (s) is not calculable within pertur- 
bative QCD (pQCD) at low energies and is ob- 
tained in a similar way by means of a dispersion 
relation. Current electroweak precision tests at 
the scale s — M§ are unfortunately limited by the 
precision of the order of 1% by which the hadronic 
piece to the running of a is known. 



Fig. [T] (from ref. ;5j) shows the contributions of 
different energy intervals of hadronic cross section 
data to the absolute value and to the squared er- 
ror for a^ adr and Aa^^Mf). We see that in the 
case of the muon anomaly, the radiative return 
program described above covers the entire en- 
ergy range of interest. The two-pion cross section 
(pion formfactor) is the by far dominating chan- 
nel, since below 1 GcV it contributes to ~ 60% to 
a'i adr . In the case of the fine structure constant 

~ 30% of the absolute value of Aa^ r (M|) are 
covered by the mass range < 4.5 GeV. If pQCD 
is applicable down to lower energies or in case 
novel theoretical techniques can be used in the 
evaluation of the dispersion integral, the impact 
of the hadronic cross section data below 4.5 GeV 
is very much enhanced. For further details con- 
cerning the evaluation of the dispersion integrals 
we refer to ref. [3] [4] [5]. 

2. THE KLOE MEASUREMENT OF 
THE PION FORMFACTOR 

KLOE has published a radiative-return anal- 
ysis of the pion formfactor, in which the ISR- 
photon is required to be emitted at small (large) 
polar angles 6 1 < 15° (0 7 > 165°) [6]. No photon 
tagging is possible in such an approach. The 
high momentum resolution of the KLOE drift 
chamber allows a separation of signal from back- 
ground with very high precision even without 
the kinematic closure of the event. A new and 
complementary analysis at large photon angles 
(50° < 9 1 < 130°) and with tagged photons 
is under study now and will allow to cover the 
threshold region, which was kinematically for- 
bidden before. Moreover, the beforementioned 
normalization to radiative muon pairs will allow 
an important cross-check of the radiator function 
approach. 

2.1. Small angle analysis 

The small angle analysis provides high statis- 
tics for ISR-photons and suppresses the relative 
amount of events, in which the photon is emitted 
from the final state pions (final state radiation, 
FSR). FSR is an irreducible background to the 
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radiative return analysis. The relative amount of 
FSR-events is well below 1% for the chosen selec- 
tion cuts and has been estimated by means of the 
Monte-Carlo-generator PHOKHARA, which uses 
the model of scalar QED (pointlike pions) as a 
model for FSR. Special attention has been given 
to events, in which simultaneously an ISR- and 
FSR-photon are emitted (NLO-FSR) since those 
events must not be considered as a background 
and are also not suppressed by the acceptance 
cuts. The relative contribution of NLO-FSR 
events is known with a precision of 0.3%. 
A total experimental error of 0.9% has been 
achieved for the e + e~ — * 7r + 7r~7 cross section 
measurement. The error consists of the individ- 
ual contributions of the selection efficiencies and 
of the precision, with which the residual back- 
ground after all selection cuts is known; further 
details can be found in ref. [6] . The pion formfac- 
tor, which can be obtained from the non-radiative 
cross section e + e~ — > ir + ir~, is extracted with 
a total error of 1.3%. This includes also the 
theory uncertainties associated with the radiator 
function and with the large-angle Bhabha cross 
section. The Bhabha cross section is needed for 
the luminosity measurement and will diminish 
with the new version of the BABAYAGA event 
generator, see ref. [7]. 

The KLOE measurement of the pion formfac- 
tor has been used to compute the contribution of 
the two-pion cross section to a^ adr in the energy 
interval 0.35 < < 0.95 GeV 2 . The value 

a™ = (388.7 ± 0.8 stat ± 3.5 syst ± 3.5 thco )xl0- 10 
covers about 60% of the total contribution. The 
KLOE result agrees within 0.5 standard devia- 
tions with values for aV' computed from the 
data sets of the experiments CMD-2 [8] and 
SND [10], which were operated in the last years 
at the VEPP-2M-collider in Novosibirsk. The 
relative difference of the mass spectra is shown 
in fig. [3] For this comparison the KLOE data 
points have been interpolated and the measured 
data points of CMD-2 and SND are used in the 
plot. We observe relatively large deviations of 




1 in the somewhat smaller mass mass region 0.37 < M. 
0.93 GeV 2 
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Figure 2. KLOE measurement of the pion formfactor 
using radiative-return events with photons emitted at 
small polar angles. 



up to some percent between KLOE and SND 
at high and low masses, while the overall agree- 
ment is better with CMD-2. Please notice that 
SND very recently has updated [5j its previous 
measurement ref. [10] . taking into account a new 
treatment of radiative corrections. The disagree- 
ment with KLOE has diminished with new data, 
but a clear trend of difference is still visible. 
The good agreement in the dispersion integral is 
partly due to a compensation effect at lower and 
higher energies. 

The KLOE measurement presented above refers 
to data taken in 2001 with a total integrated lumi- 
nosity of ~ 140 pb _1 . KLOE is now performing 
an analysis using 2002 data (~ 240 pb _1 ) [11], for 
which a total systematic error (experimental and 
theoretical) < 1% is expected. The acceptance 
cuts of this analysis will be unchanged; the im- 
provement will be due to better and more stable 
running conditions as well as due to modification 
in the online and offline environment, which will 
result in lower systematic errors associated to the 
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Figure 3. Relative difference of the pion formfactor 
measurements from the experiments CMD-2 (trian- 
gles) and SND (circles), relative to the KLOE mea- 
surement. The KLOE data points have been interpo- 
lated and the statistical (light grey) and the system- 
atic (dark grey) error bands are shown in the plot. 



trigger and background-filter efficiencies. Con- 
cerning theory, a new version of BABAYAGA is 
available [7] , which allows to reduce the luminos- 
ity error by about a factor 2. The final goal of the 
analysis is a measurement of R, which requires 
that e + e~ — ► fi + fi~j events are selected with 
equally high precision as 7r + 7r _ 7 events . 

2.2. Large angle analysis 

The analysis described above, in which the 
ISR-photon is emitted at small polar angles, 
does not allow to cover the threshold region 
< 0.35GeV 2 , since in this kinematical region 
the two pions are emitted essentially back-to-back 
to the ISR-photon and hence cannot be detected 
simultaneously in the fiducial volume defined for 
the pion tracks 50° < 6^ < 130° [12]. In order to 
measure the pion formfactor at threshold, KLOE 
is now performing a complementary analysis, in 



which the ISR-photon is tagged at large polar 
angles 50° < 9 7 < 130°. Due to the 1/s 2 de- 
pendence in the dispersion integral for ajj adr , the 
low mass region of the two-pion cross section is 
actually giving a ~ 20% contribution to the total 
integral and hence an improved determination of 
the cross section at threshold is needed. 
At large photon angles background from (f> — > 
ir + ir~ir° is huge and dedicated selection cuts, like 
a cut on the angle between the missing momen- 
tum and the tagged photon direction, as well as a 
kinematic fit in the background hypothesis with 
a cut on Xtttvk are needed to suppress this contri- 
bution. Moreover, irreducible background from 
events with the same 7r + 7r~7 final state is not 
negligible anymore at large photon angles. This 
background category, which has to be subtracted 
relying on Monte-Carlo prediction, is given by 
FSR-events and by the (j> radiative decay into 
the scalar / (980) with / (980) -> tt+tt-. A 
possible model dependence of FSR (the model 
of scalar QED is used in PHOKHARA) and of 
the description of the scalar f (980)7 amplitude, 
can be tested by means of the forward-backward 
asymmetry. In a recent KLOE publication [13] 
good agreement between data and simulation has 
been found for the forward-backward asymme- 
try, setting upper limits for the systematic errors 
associated with these model uncertainties. 

The main limitation for the measurement of 
the pion formfactor at threshold will arise from 
the background channels <f> — ► 7r + 7r~7r° and 
<f> — ► /o(980)7 — > 7r + 7r _ 7. In order to further 
reduce the systematic errors associated to these 
channels, the DA<I>NE collider has taken data off- 
resonance at a center-of-mass energy of y/s = 1.00 
GeV in its last KLOE run (December 2005 to 
March 2006, 250pb _1 integrated luminosity). 
The off-peak analysis will allow a considerably 
improved determination of the threshold region. 
Moreover, together with the data taken on-peak 
it will be possible to study the interference of the 
/o(980) amplitude with FSR. 
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3. THE RADIATIVE RETURN AT PEP- 
II WITH THE BABAR DETECTOR 

The BaBar experiment has previously 
published results of the exclusive 7r + 7r~7r , 
7r + 7r~7r + 7r~, K + K~tt + it~ and K + K~K + K~ 
final states with better coverage than all previous 
experiments and comparable or better preci- 
sion, using 89 fb _1 of data. New results are 
presented here for final states with 6 hadrons 
and for the timelike proton-antiproton formfac- 
tor. The new results contain a larger data set 
of ~ 240pb~ 1 . The mass range of interest in all 
these measurements Mh a( j r < 4.5 GeV is rather 
distant from the center-of-mass energy of PEP-II 
(^/s = m X ( 4S ) = 10.58 GeV), requiring hence 
a high-energetic ISR-photon with an energy in 
the center-of-mass system of E* > 3 GeV. Such 
an event signature separates ISR-events from 
background from the T(4s)-resonance (B-decays 
and their consecutive decay products), since pho- 
tons originating from these events are in general 
much less-energetic. The event signature at the 
B-factory also implies the necessity for photon 
tagging in order to have high acceptance for 
the hadronic system. The mass resolution is 
improved by performing a kinematic fit in the 
signal hypothesis, requiring constraints on four- 
momentum and (in case of neutral pions) of the 
7r°-mass. A cut on the \ 2 of the kinematic fit 
is the main tool for background subtraction and 
the shape of the ^-distribution is used with a 
sideband-method to obtain the residual back- 
ground after all cuts. 

In this paper we report not only on the cross 
section measurements, but we also present stud- 
ies of the internal structures. For all channels 
discussed here the J/ip signals were used to 
extract the associated branching fractions and 
in some cases this has been done also for the 
ip{2S) signals. BaBar has published 10 J/ip and 
3 ip(2S) branching ratios; in 4 cases there has 
been no previous measurements, in 5 cases the 
measurements are better than the current world 
average. In the following we give a brief overview 
of the individual channels published so far by 
BaBar with special emphasis on the more recent 
analyses of 6 Hadrons and e + e~ — > pp. 
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Figure 4. BaBar measurement (filled circles) of the 
cross section e + e~ — > 7r + 7r _ 7r°, compared to previ- 
ous measurements from SND (open circles) and DM2 
(open triangles). 



3.1. e+e" -> 3 Pions 

The 7T + 7r~7r° mass spectrum has been mea- 
sured from 1.05 GeV up to the J/ip mass region 
with a systematic error of ~ 5% below 2.5 GeV 
and up to ~ 20% at higher masses [14]. The 
spectrum is dominated by the lo, <p and J/ip res- 
onances. The BaBar measurement could improve 
significantly on the world's knowledge of the ex- 
cited lo states. Fig.|4]shows the BaBar data points 
together with previous measurements in the en- 
ergy range relevant for the exited states to' and 
lo". A clear disagreement with DM2 [TS] is seen 
above 1.4 GeV. The BaBar spectrum has been fit- 
ted up to 1.8 GeV and the following results for the 
masses and widthes of the lo' and lo" have been 
foundl M{lo') = (1350 ± 20 ± 20) MeV, T(lo') = 
(450 ± 70 ± 70) MeV, M(lo") = (1660 ± 10 ± 2) 
MeV, T{lo") = (230 ± 30 ± 20) MeV. 

3.2. e+e — > 4 Hadrons 

The 7r + 7r~7r + 7r~, K + K~tt + t:~ and 
K + K~ K + K~ exclusive final states have been 
measured from theshold up to 4.5 GeV with 
systematic errors of 5%, 15% and 25%, respec- 
tively [IB]. Fig. [5] shows the mass distribu- 
tion of the events satisfying the 7r + 7r _ 7r + 7r~ and 
K + K~tt + tt~ kinematics as well as K/tt identifi- 

2 in this fit M w and the phases between to, <j> with respect 
to io', lo" are fixed. 
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Figure 5. The energy dependence of the 
e + e~ — > 7r + 7r _ 7r + 7r _ (upper plot) and the e + e~ — » 
A' + if _ 7r + 7r^ (lower plot) cross section obtained by 
BaBar (filled circles) by radiative return in compari- 
son with previous data. 



cation constraints. We identify that the radiative 
return technique at BaBar allows to cover a 
wide energy range in one single experiment. The 
K + K~ K + K~ measurement is the first measure- 
ment ever. Background is relatively low for all 
channels under study (e.g. few percent at 1.5 
GeV for 7r + 7r - 7r + 7r~) and is dominated by ISR- 
events of higher multiplicities and of continuum 
non-ISR events at higher masses. 
The 7T + 7r~7r + 7r~ final state is dominated by 
the two-body ai(1260)7r intermediate state; the 
K + K ~Tr + ir~ final state shows no significant 
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Figure 6. The energy dependence of the e + e~ — > 
3(7r + 7r~) (upper plot) and the e + e" — » 2(iy + ty~)2tt 
(lower plot) cross section obtained by BaBar (filled 
circles) by radiative return in comparison with previ- 
ous data. 



two-body states, but rich three-body struc- 
ture, including K* (890)Kir, <j>Trir, pKK and 
K*(U30)Kir. 

3.3. e + e _ — > 6 Hadrons 

The 6-Hadrons final state has been measured 
in the exclusive channels 3(7r + 7r~), 2(tt + tt^)2tt° 
and K + K~2(ir + ir~) [17]. The cross section in the 
last case has never been measured before; the pre- 
cision in the first two cases is ~ 20%, which is a 
large improvement with respect to existing data. 
Again, the entire energy range from threshold up 
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Figure 7. BaBar measurement of the ratio of the 
electric and magnetic formfactor describing the cross 
section e + e~ — > pp. Filled circles depict BaBar data, 
the curve is a fit result. Open circles show data from 
the experiment PS170. 



to 4.5 GeV is measured in one single experiment. 
The distributions for the final states 3(7r + 7r~) and 
2(7T + 7r-)27r° are shown in fig. [6l A clear dip 
is visible at about 1.9 GeV in both modes. A 
similar feature was already seen by FOCUS [TH] 
in the diffractive photoproduction of six charged 
pions. The spectra are fitted using the sum of 
a Breit-Wigner resonance function and a Jacob- 
Slansky continuum shape. For the 3(7r + 7r~) 
(2(7T+7r-)27r ) mode, we obtain values of 1880±30 
MeV (1860 ± 20 MeV) for the resonance peak, 
130 ± 30 MeV (160 ± 20 MeV) for the resonance 
width and 21° ± 14° (-3° ± 15°) for the phase 
shift between the resonance and continuum. The 
width values differ significantly from the FOCUS 
result of 29 ± 14 MeV. 

Also the substructures of the individual modes 
have been investigated. There is surprisinlgly lit- 
tle substructure in the 3(7r + 7r~) channel. The 
spectrum can be described by a simulation with 
one p° and four pions distributed according to 
phase space. The 2(7r + 7r~)27r° final state shows 
a similar cross section but a much more com- 
plex internal structure. BaBar observes signals 



for p° , p^, ui and r/, and a substantial con- 
tribution from the two-body uir] intermediate 
state, which appears to be resonant. Additional 
modes containing more tt° are under study. Also 
the K + K ~2(7r + 7r~) channel shows an interesting 
substructure with a rather weak <p contribution, 
but a strong if* (890). 

3.4. e + e~ — > pp 

BaBar has also performed a measurement of 
the e + e~ — > pp cross section [20] . The experimen- 
tal challenge in this case is to fight the KK, irir 
and pp, two-body background, which has a cross 
section approximately a factor 10 — 100 higher. 
Apart from the good momentum resolution of the 
detector, BaBar has excellent particle identifaca- 
tion capabilities, which allow a clean selection of 
protons; ~ 4000 pp events are selected from a 
data sample of ~ 240pb~ . The timelike proton- 
antiproton form factor is parametrized by an elec- 
tric and a magnetic formfactor Ge, Gm- 
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where the factor C accounts for the Coulomb in- 
teraction of the final state particles. The proton 
helicity angle 8 P in the pp rest frame can be used 
to separate the \Ge\ 2 and | G a/ | 2 terms. Their 
respective variations are approximately ~ sin 2 8 p 
and ~ (1 + cos 2 (9 p ). By fitting the cos8 p dis- 
tribution to a sum of the two terms, the ratio 
\Ge/Gm\ can be extracted. This is done sepa- 
rately in six bins of M p p. The results are shown 
in fig. [71 and disagree significantly with previous 
measurements from LEAR 21} close to threshold. 
At larger values of M pp the BaBar measurement 
finds \G E /G M \ ~ 1. 

In order to compare the cross section measure- 
ment with previous data (e + e~ and pp experi- 
ments), the effective form factor is introduced: 

G = J\G E \ 2 + 2ml/s\G M ?- The BaBar mea- 
surement of G is in good agreement with existing 
results, as can be seen in fig. [5J The structure 
of the formfactor is rather complicated; we make 
the following observations: (i) BaBar confirms an 
increase of G towards threshold as seen before 
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Figure 8. The e + e~ — > pp cross section measured by 
BaBar (filled circles) in comparison with data from 
other e + e~ colliders. 



by other experiments; (ii) two sharp drops of the 
spectrum at M pP = 2.25 and 3.0 GeV are ob- 
served; (iii) data at large values M pP > 3 GeV is 
in good agreement with the prediction from per- 
turbative QCD. 



4. CONCLUSIONS 

At the particle factories DA$NE and PEP-II 
the radiative return method allows precision mea- 
surements of the hadronic cross section, which are 
of utmost importance for an improved determi- 
nation of the hadronic contribution of the muon 
anomaly a^ adr and of the running fine structure 
constant a(M§). At DA$NE the KLOE exper- 
iment has measured the pion formfactor with a 
precision of 1.3%. An update of the analysis will 
lead to a further reduction of the systematic er- 
ror and a normalization to radiative muon pairs 
is foreseeen. Moreover, KLOE has collected data 
with the collider running off-resonance, which will 
allow an improved determination of the threshold 
mass region. 

At the B-factory PEP-II a program to measure 
all exclusive hadronic channels from theshold up 
to 4.5 GeV is underway. In addition to the pub- 
lished results presented in this paper, there are 
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ongoing analyses for further final states, includ- 
ing the important two-pion- channel. 
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